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Abstract. We analyze the quantum dynamics of a non-relativistic particle moving 
in a bounded domain of physical space, when the boundary conditions are rapidly 
changed. In general, this yields new boundary conditions, via a dynamical composition 
law that is a very simple instance of superposition of different topologies. In all 
cases unitarity is preserved and the quick change of boundary conditions docs not 
introduce any decoherence in the system. Among the emerging boundary conditions, 
the Dirichlct case (vanishing wave function at the boundary) plays the role of an 
attractor. Possible experimental implementations with superconducting quantum 
interference devices are explored and analyzed. 
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1. Introduction 



The dynamics of a quantum particle in a bounded domain of physical space is a 
paradigmatic problem in quantum mechanics. From the mathematical point of view, 
one must make sure that physical observables be properly defined in terms of self- 
adjoint operators [1, 2]. This translates into suitable choices of boundary conditions, 
those most commonly used in physics being Dirichlet (vanishing wave function at the 
boundary), Neumann (vanishing normal derivative), and periodic ones. From a physical 
perspective the problem has a plethora of applications, ranging from atoms in cavities 
[3] to ions and atoms in magnetic traps [4], to superconducting quantum interference 
devices (SQUID) [5]. The role and importance of boundary conditions has been recently 
stressed in an interesting article [6], where varying boundary conditions are viewed as 
a model of spacetime topology change. Notable applications arise in the context of the 
Casimir effect and its dynamical version, giving rise to photon generation in a microwave 
cavity with time-dependent boundary conditions [7]. In this Article we shall analyze the 
effect of a rapid alternating change of boundary conditions and show that the resulting 
dynamics yields novel boundary conditions. We shall make use of a general result [8], 
that characterizes all possible unitarity-preserving boundary conditions. 



2. Trotter formula for alternating boundary conditions 

Let us consider a spinless particle of mass m in a cavity Q with a regular boundary dQ. 
We focus on the dynamics that arises when the cavity undergoes a rapid alternating 
change of boundary conditions with a time period t/N. Such an evolution is described 
by the following unitary operator 

(^Q-itTu/N^-itTy/N^ ^^-itTu/N^-itTv/N-^ _ _ _ (^^-itTu /N ^-itTy /N^^ 

' V ' 

A'' times 

where Tu = /2m = —Au/2'm is the kinetic energy operator, with Ajj a self-adjoint 
extension of the Laplacian A, with given boundary conditions specified by the unitary 
operator U acting on the boundary Hilbert space L^{dQ) [8]. Notice that in general A[/ 
and Ay do not commute. Physically, this corresponds to a rapid switching between two 
different boundary conditions, in the limit of frequent switching. 

Evolutions of this kind are familiar in the context of quantum chaos [9, 10] and in 
connection with the quantum Zeno effect [11, 12] (when the von Neumann measurements 
[1] of the most familiar formulation d la Misra and Sudarshan [13] are replaced by 
frequent unitary pulses, as in the experiment by Itano et al [14].) The analogy between 
these apparently unrelated dynamics has been explored by several authors [15]. 

The relevant question is to show whether in the N ^ oo limit (when the time 
interval between the switches goes to zero, the number of switches goes to infinite, while 
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the total time 2t is kept constant), the evolution is given by 

in terms of a Hamiltonian Ti^, with some boundary conditions W . 
We will prove that this is the case, and in fact 

W = Ui.V = V ^U, (3) 

where * is a commutative and associative product on the boundary unitary operators 
defined by 

U,V:=4^^^H)±^\. (4) 



2 

Here €. is the Cayley transform that maps Hermitian into unitary matrices: 

^(^) - ^-(^) - (^) 

Notice that the Cayley transform is not onto. Its range is the subset of unitary matrices 
whose eigenvalues are different from —1. Thus, strictly speaking, Eq. (4) has a meaning 
only for U and V in this subset. We will show what is the action of -k on all boundary 
unitaries and prove that the eigenspaces with eigenvalues —1 are absorbing for the 
product. In particular, 

(_/) ^v = -I, (6) 

for any unitary V . 

The solution makes use of some mathematical results on product formulae. In a 
seminal paper [16] Trotter proved that if A and B are self-adjoint and C = A + B is 
self-adjoint on the intersection of their domains, D{C) = D{A)r\D{B), then the formula 

holds for — > oo. This is the famous Trotter's product formula. 

Unfortunately, this formula cannot be applied to our case. Indeed, on the 
intersection of the domains of the Laplacians with different boundary conditions, 
D = D{Tij) n D(Tv), the Laplacian is not self-adjoint since the domain is too small, 
being defined by too many constraints (those of U and those of V). Even more, it 
admits many self-adjoint extensions, so that the meaning of Tw in (2) is unclear: which 
boundary conditions W, if any, are obtained in the limit? 

The answer is obtained by considering the quadratic forms associated to the 
operators (i.e. their expectation values), instead of the operators themselves. Notice, in 
fact, that the domain D{tu) of the quadratic form of the kinetic energy, tjj^ip) = {jp, Tjjip) 
(see below), is larger than the operator domain D{Tu), and there is a unique self- 
adjoint kinetic operator T\y associated to the quadratic form {tu + tv)/2 on the (dense) 
intersection D = D{tu) D D{tv)- The operator Ty/ defined on D{Tv/) d D is called the 
form sum of Tjj and Ty and is denoted by Tw = (Tu+Tv)/2. 

This idea, introduced by Kato [17], was elaborated by Lapidus [19, 18] who found 
the ultimate version of Trotter's product formula (7): If A and B are self-adjoint and 
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bounded below, and the intersection of their form domains D{a) fl D{b) is dense, then 
Trotter's formula (7) holds when — )■ oo, with 

C = A+B, (8) 

the form sum of A and B. 

As a technical remark notice that, as a consequence of the weakening of the 
hypotheses, the convergence of the product formula when the operator sum is not self- 
adjoint is in a weaker topology than in Trotter's case. More precisely, if C = A + B 
is self-adjoint, then the convergence is pointwise (in fact locally uniform) in t in the 
strong operator topology. On the other hand, when C = A+B the limit is proved only 
on the average in t. Whether this convergence result is physically satisfactory or not, 
and whether it can be made stronger is a long-standing problem. For an interesting 
discussion see [20]. 



3. Boundary conditions for a free particle on an interval 

In order to simplify our analysis and focus on concrete physical problems we shall restrict 
our attention to the problem of a particle moving in the interval fl = [0, 1] of the real line 
and set for convenience m = 1/2. Apart from being the simplest mathematical case, the 
one-dimensional problem is also interesting in its own. Indeed, it can be implemented in 
a SQUID with a tunable junction, obtained by replacing the junction with an additional 
flux loop [21, 22, 23]. On such devices different boundary conditions can be implemented 
by tuning the magnetic flux across the additional loop. More generally, a suitable 
combinations of tunable SQUIDs can be an experimental realization of the continuous 
interpolation among different topologies of the model introduced in [6]. 

The whole family of self-adjoint extensions of the Laplacian is in one to one 
correspondence with the possible boundary conditions given by the 2x2 unitary matrices 
of [/ G W(2) [8]. They are given by the Laplacian 

d^ 

acting on the domain D(Tu) of functions belonging to H^{0, 1), the Sobolev space of 
square integrable functions ip with square integrable second derivative ip", and satisfying 
the boundary conditions 

i{I + U)ip' = {I - U)^, U eU{2), (10) 



where 





^ - " . V := r . (11) 



Physically, U = — / corresponds to Dirichlet boundary conditions, 

^(0) = V^(l) = 0, (12) 
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while Neumann boundary conditions, 

^'(0) = ^'(1) = 0, (13) 



are given hj U = I. The latter is a particular case of Robin's boundary conditions 
given by U 



ry (y 

i/.'(0) = -taii-^i(0), i/.'(l) = tim-^i(l). (14) 



I (15) 

corresponds to Dirichlet at the left and Robin at the right, 

V;(0) = 0, ^'(l)=tan-V^(l). (16) 

For nondiagonal U the boundary conditions at the two ends mix. For example, pseudo- 
periodic boundary conditions, 

V^(l)=ei"^(0), ^A'(l)=e'V(0), (17) 

are obtained by 



e" 




U = \ = COS a ar, + sin a ay. 



In general, if the unitary U has no —1 eigenvalues, the wavefunction can assume 
any values at the ends. Only the boundary values of its derivative are constrained is 
some way. These unitaries, corresponding to free ends, will be called regular. Other case 
are not regular. For instance, one eigenvalue equal to —1, as in (15) or (18), corresponds 
to one constraint on the values of at the end, as in (16) or (17), respectively. Finally, 
two —1 eigenvalues, i.e. U = —I, correspond to two constraints on the wavefunction 
boundary values (12). 



4. Quadratic forms 

Let T[/ be the kinetic energy operator. For any i/j G D{Tjj) C i^^(0, 1), an integration 
by parts gives 

tu{il^) = {tp.Tuij) = - [ tP{x)tP"{x)dx 

Jo 

= j \i)'{x)\^dx-i}{i)i)'{i) + ^{Q)i)\Q) 

= ||^'f -(^1^'). (19) 

We now try to rewrite the boundary form in a more convenient way by making use of 
the boundary conditions of Tjj. In particular we want to trade the boundary values of 
the derivative for the boundary values of the function, in order to obtain a quadratic 
form 

tu{i^) = Wf-Tu{^). (20) 
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with Tu{if) a quadratic form of the boundary vector ip, given in (11). In this way, 
the quadratic form (20) is defined on the Sobolev space H^{0,1) of square integrable 
functions with square integrable first derivative, at variance with the operator Tjj which 
required also the square integrability of the second derivative. Notice indeed that, since 
■ip' is in L^(0, 1), then il^lx) = 'ilj'{y)dy is a continuous function and its boundary 
values (p are well defined. On the other hand, the boundary values of its derivative ip' 
would make no sense. 

We distinguish among three possibilities according to the number of eigenvalues of 
U equal to —1. Let 

U = u,\0{^\ + u,\^^){^^l (21) 

with = \u2\ = 1, and (^|^^) = 0, its spectral decomposition (here {^\ri) = 
^iVi + ^2V2)- We have: 

(i) If U is regular, i.e. ui^2 7^—1, then (J + U) is invertible, and the boundary values 
of the derivative can be expressed in terms of the boundary values of the function 

^' = Kuif. (22) 

where Ku = €~^{U) is the inverse Cayley transform of U, defined in (5). Therefore, 
the boundary form F^/ is the expectation value of Ku, namely 

Tuiv) = (vlKuif), (23) 
and D{tu) = H^{0, 1), with free ends (no constraints on the boundary values ip). 

(ii) If —1 is a nondegenerate eigenvalue of U, that is ui = —1 and U2 7^ —1, then U is 
not in the range of the Cayley transform (5), but from (10) we get (^Iv^) = 0, and 

n^') = -i^nv^)- (24) 

Therefore, 

m = 0, ru{^) = i\^\{e\^)\\ (25) 

I + U2 

(iii) Finally, if Mi = U2 = —1, then U = —I, and 

ip = 0, T_j{ip) = 0. (26) 

5. Composition law of boundary conditions 

We now evaluate the limit of the alternating dynamics (2). According to (8) the product 
formula (2) holds with the form sum 

Tw = l [Tu+Tv] . (27) 

Thus, the evaluation of the emergent boundary condition W in (2) requires the 
computation of the sum 

tw = ^ {tu + tv) , (28) 
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and its domain 

D{tw) = D{tu)nDitv). (29) 
Again, we distinguish various cases according to the number of eigenvalues —1: 

(i) In the regular case, when both U and V have no eigenvalues equal to —1, we get 
from (22) that = {Kjj + Ky)/2, with no constraints on the wave-function 
boundary values if. Therefore, one can write W = U -kV with U -kV given by (4). 
Explicitly, we get 

/ _ 1 fl^ + lull) 

w = u.v= (30) 

2 \I+U I+V) 

(ii) If —1 is a nondegenerate eigenvalue of f/, and V is regular, then D{tw) = D{tu), 
with the only constraint = 0. Therefore, the boundary forms Tu and Ty are 
nonzero and add up only on the orthogonal subspace, spanned by It is easy to 
see that 



W = U^V = -\0{^\+W2\e){^\ (31) 



1-i 



W2 = 7 (32) 

1 ( l-«2 I /C±\I-V, ^ 



with 



1 / 1-U2 I /C±\ J-V C±\ 

2 \1+U2 ~'~ I I+V'^ I 

(iii) If —1 is a nondegenerate eigenvalue of both U and V, that is mi = wi = — 1 and 
W2,'i^2 7^ —1, then there are two possibilities 

(a) If the eigenvectors of U and V belonging to —1 are parallel, that is U commutes 
with V, then D(tw) = D(tu) = D(tv)- Thus, the only constraint is (^1^9) = 
and W has the previous form (31), where (32) particularizes into 

1 _ 1 I 1""2 _|_ l"<^2 \ 
2 \l+U2 ^ 1+V2 

W2 = 7 ^. (33) 



H 


41 




—U2 




-V2 \ 








+U2 




t-f2 / 



(b) If the eigenvectors ^ of f/ and 77 of belonging to —1 are not parallel, then 
they span the whole space. The constraints (^Iv^) = and {Ti\ip) = imply 
Dirichlet's boundary conditions = 0, so that D{t]v) = D{t-j) and 

W = Ui.V = -I. (34) 

(iv) Finally, in the case U = —I (or V = —I) then D{tw) = D{t^i), so that ip = and 

= (-/)* t/ = t/ *(-/) = -/. (35) 



6. Prediction: a SQUID circuit 

Summarizing, the limit boundary condition is W = U -kV, where -k is given in terms of 
the Cayley transform (4) for regular U and V (free ends (f), while all constraints on the 
wave-function boundary values are conserved by the product and inherited by W. 



Composing boundary conditions 



8 



Therefore, in one dimension, if U and W have independent constraints on (p, then 
Trotter yields Dirichlet boundary conditions, (/? = 0. A particularly interesting case 
where this happens is for pseudo-periodic boundary conditions, where U and V have 
the form (18) with two different phases ai ^ a2- In that case the boundary conditions 
read 

^(0) = e^"^^(l) = e'^^-^il), (36) 

and imply that ^/'(O) = ^(1) = and W = —I. This case can be experimentally 
implemented by means of a SQUID circuit where the properties of the Josephson 
junction are pulsed to mimic the Trotter evolution described above. The result 
(somewhat counterintuitive) will be a blockage of the electrical current through the 
circuit. 

This situation is very different from the seemingly similar problem where the 
properties of the Josephson junction is not modified, but a magnetic flux across the 
SQUID is pulsed in a similar way from ai to a2- In this case the resulting evolution is 
very different because the boundary conditions are always the same but the Hamiltonians 
are not the Laplacian but have magnetic couplings 



In this case the Hamiltonian form sum is 

dx 



^^3 = ^ = + "3 , (3J 



where 



as = (39) 



7. Conclusions 



We have analyzed the quantum dynamics of a non-relativistic particle moving in a 
bounded domain of physical space, when the boundary conditions are rapidly changed. 
We have seen that the resulting boundary conditions is obtained though a dynamical 
composition law that can be viewed as a simple instance of superposition of different 
topologies. Having in mind quantum gravity and the associated fluctuating topologies, 
one may speculate that according to our composition rule, the Dirichelet boundary 
conditions (that act as an attractor) play the same role as the classical paths in Feynman 
path integral, i.e. they completely dominate in the classical transition. 
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